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Abstract Hypothesis of quark binding through condensation of gluons inside hadrons is formulated 
in the context of a renormalization group procedure for effective particles (RGPEP) in the light- 
front (LF) Hamiltonian approach to QCD. At the momentum scales of relative motion of hadronic 
constituents that are comparable with Aqcd, the hypothetical boost-invariant constituent dynamics 
is identified using gauge symmetry. The resulting picture of mesons and baryons closely resembles 
constituent quark models with harmonic oscillator potentials, shares some features of AdS/QCD, and 
can be systematically studied using RGPEP in QCD. 
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1 Canonical approach to QCD 

The canonical plan for solving QCD is to start from the corresponding classical gauge-invariant La- 
grangian density, C can , derive from it a candidate for a Hamiltonian density, T-L C an, and evaluate the 
Hamiltonian H can — f d 3 x T-L ca n , in which the quark and gluon fields are quantized. Our discussion of 
the quantum theory begins with the standard form of Hamiltonian dynamics, called the instant form 
(IF) , but our goal is to address the front form (FF) of the theory pQ . 

The canonical operator H can acts in the Fock space of states \tp) that are built from a vacuum by 
the operators that create quarks and gluons. Then, a solution to the eigenvalue problem 

H can \ 1 p)=E\^), (1) 

with energy of the form E = ^JM\ + P 2 , where M h is a mass, should represent a hadron h moving 
with momentum P. Time evolution of all states would be described using operators of the form 
U{t2,t\) = e~ lHcan ( t2 ~ tl \ etc. The ultimate goal of such plan for QCD is to achieve a quality of 
the wave function representation of hadrons that matches, and in the future hopefully even exceeds 
the quality of today's QED representation of atoms and their chemistry. 



2 Time-honored problems of canonical approach 

The canonical approach encounters divergence problems that are known for a long time to require a 
logical resolution [5] . A hint of the resolution was discovered also a long time ago [5] through the concept 
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of a renormalization group (RG) procedure [1]. The idea is that one can calculate the counterterms 
required by the adopted regularization. This is done in a sequence of steps that eventually produce a 
manageable effective theory. To begin with, one has to replace H can [5J E] with a regulated operator, 
Hcam where A stands for an extreme cutoff parameter. The counterterms, CT A , that need to be 
included in the initial Hamiltonian, H = H A an + CT A , are obtained from the condition that the 
matrix elements of the desired effective Hamiltonian do not depend on the regularization in H . The 
calculation of counterterms is based on evaluation of a whole family of effective Hamiltonians, H\ 
labeled by the RG parameter A, which can be chosen to have the dimension of momentum. This 
parameter plays the role of a sliding momentum cutoff. 

The question that so far has no satisfactory resolution in canonical approaches to QCD is: What 
precisely is the operator H\ in which A is so small that the eigenvalue problem for H\ can be solved 
using computers? One well-known reason for the difficulty is asymptotic freedom [7] [S] . It implies that 
the effective coupling constant g\ in H\ increases according to the rule g\ ~ 1/ \n{\ / Aqcd) when A 
decreases (for the case of H\ in LF QCD, see [9]), and becomes too large too soon for a perturbative 
procedure to produce H\ with a sufficiently small A for a reliable computation of quark and gluon 
wave functions of hadrons. 

The chief difficulty of the IF of dynamics, however, is that we do not know the ground state of 
QCD. The canonical IF interaction Hamiltonians are able to create virtual particles from empty space 
and as time flows they generate this way an infinitely complex state that so far nobody can describe 
[2]. The ground state, called vacuum, should be invariant with respect to the Lorentz transformations. 
But such transformations change momenta by arbitrary amounts and a theory with a finite cutoff 
on momentum cannot have a Lorentz invariant ground state. Nevertheless, a covariant perturbation 
theory in QCD allows for introduction of non-trivial parameters that can be interpretated as vacuum 
expectation values of operators. These parameters are related to the spectrum of hadrons through 
dispersion relations. This is how the quark and gluon vacuum condensates are introduced in the QCD 
sum rules [TU1 E] • 

In contrast to the case of the IF of dynamics, it is well-known that the FF of Hamiltonian dynamics 
necessarily leads to a trivial vacuum state, apparently allowing one to make progress without getting 
stuck in the hard vacuum problem right away. This fact motivated efforts to set up a similarity renor- 
malization group (SRC) procedure for Hamiltonians [12l [13] and apply it to LF QCD [14]. The effects 
that in the IF are associated with the unknown state of the vacuum, were suggested in Ref. [14] to be 
contained in the new terms in H, and thus also in H\, that a precise SRG procedure could identify. 
Through later efforts, summarized in Ref. [15] . it was found that there exists an additional possibility 
for the vacuum-like terms to emerge in LF QCD in a systematic calculation. The additional possibility 
is the subject of this article. 

3 Renormalization group procedure for effective particles 

In the renormalization group procedure for effective particles (RGPEP) that is used in Ref. [13] to 
suggest a new locus for vacuum-like condensate effects in LF QCD, the key element of reasoning is a 
unitary connection between creation and annihilation operators for bare, canonical quarks and gluons, 
and the operators for effective particles, 



The parameter s = 1/A has interpretation of the size of the effective particles with respect to strong 
interactions. The parameter A describes the momentum width of vertex form factors in interaction 
terms, and s, as its inverse, corresponds to the non-local interaction vertex range in space. Thus, the 
subscript refers to the point-like particles, with only local interactions, while s > implies non-local 
interactions [16] . Correspondingly, the effective quantum fields in space are constructed according to 
the well-known Fourier superposition rule [15] 



a. 



'8 



U s a Ul . 



(2) 




(3) 



As a result, ipo (x) corresponds to canonical fields, and tp s (x) to fields of effective quanta of size s. 
By writing ip s (x) = iJj(x,s), one can realize that the scale-dependent effective theories form together 
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a single 5-dimensional theory, the 5th dimension being the size of effective particles. It is natural to 
expect that this size is dynamically limited from above by sqcd ~ 1/Aqcd an d that sqcd corresponds 
to the depth of bulk penetration by matter fields in AdS/QCD models. We shall comment on this issue 
near the end of the article. 

The transformation U s in Eq. ^ is constructed in such a way (see [TS] and references therein for 
computational details) that the Hamiltonian is not changed, H s (a s ) = H(clq). But the Hamiltonians 
differ in their structure: H s (a s ) is a combination of products of operators a s with coefficients c s that are 
different from the coefficients Co of corresponding products of operators a,Q in H(ao). RGPEP provides 
differential (or algebraic) equations that produce expressions for the coefficients c s . From the equality 

H s (a ) = U$H(a )U s , (4) 

and the condition Uq = 1, one obtains 

^H s (a ) = [G s ,H s (a )], (5) 

with the generator 

G s = -Ut^U s (6) 

and initial condition Ho(ao) = H(ao) = H^ an + CT A . The use of s 4 is due to dimensional reasons since 
G s is designed to have dimension of mass to power 4. The non-perturbative, boost invariant generator 
G s for RGPEP is given in Eq. (C.l) of Ref. [15]. It is expressed in terms of the operators ao in the 
form of a commutator, 

G. = [H free) H+], (7) 

where Hf ree denotes the part of H(ao) that involves only products of the form ajao (one particle oper- 
ators) . is the remaining part of the Hamiltonian. The superscript + indicates that the coefficients 
c s in each and every term in Hf are multiplied by the square of the total + momentum of the particles 
that participate in the interaction described by a given term (see Eqs. (C.2) and (C.3) in [15 ). As 
a result, when s increases, one obtains from Eq. ^ the effective interactions that are increasingly 
tempered by the vertex form factors that limit the changes of invariant masses of interacting effective 
particles. For example, solving Eq. ([5| in a lowest order in powers of the interaction strength for any 
interaction term V, one obtains 

(m(s)\V s (a s )\n(s)) = e -* 4 (A<? m -A<?„) a ( m (s)\V (a s )\n(s)) , (8) 

where \m(s)) and \n(s)) denote arbitrary states in the Fock space built using operators a s and M.i m and 
Min denote the total invariant masses of only these subsets of effective particles in the corresponding 
states that are directly involved in the matrix element of the interaction V. The RGPEP vertex form 
factors suggest that the effective LF Fock space description of hadronic states may actually converge. 
This expectation needs to be verified by explicit calculations, which is a serious challenge. 

Apparently similar to Eq. ([5| , beautiful flow equations have been developed by Wegner [T7J [TS] ■ 
The five main ways RGPEP differs are: (1) The transformation of a Hamiltonian in RGPEP is limited 
to transformations of creation and annihilation operators, which is a narrower class of transformations 
than rotating matrices, since matrices that result from evaluating matrix elements of linear combina- 
tions of products of rotated creation and annihilation operators form only a subset in the set of all 
matrices of interest in quantum mechanics; (2) Eq. ^ does not require specification of the full diagonal 
matrix elements of the evolving Hamiltonian, as Wegner's equation does; (3) The Hamiltonian as an 
operator is not altered at all by Eq. ([s]) since the coefficients c s evolve in a way that is compensated by 
the evolution of operators a s , but the RGPEP derivation of counterterms alters the initial condition 
at s — 0, a feature that Wegner's equation does not include, since it treats the initial Hamiltonian as 
given (the calculation of H is one of two generic goals of RGPEP, as it is in the SRG jT2J [T3] , the other 
goal being the evaluation of corresponding effective Hamiltonians with s sufficiently large so that they 
can be used in numerical computations, and the latter goal is shared with Wegner's flows); (4) The 
RGPEP operator calculus renders coefficients c s in the Hamiltonians H s (a s ) that can be applied to 
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Fig. 1 Visualization of the idea that a three-quark configuration in a proton evolves with the RGPEP scale 
parameter s: (a) the picture at s much smaller than the scale s c that corresponds to the constituent model, 
(b) s somewhat smaller than s c , and (c) s comparable with s c . The key consequence of the RGPEP is that the 
slow effective quarks must be large. For s ~ 1/Aqcd, they are as large as the proton itself. In the case (c), 
they overlap a lot and balance color to zero inside a large part of the proton volume. 



arbitrary states in the Fock space, instead of only to a specified set of states that is used in defining the 
Wegner matrix equations; (5) The generator G s in Eq. Q is constructed to preserve the 7 kinematical 
symmetries of the FF of Hamiltonian dynamics. The feature (5) is essential for application of RGPEP 
to QCD since it allows one to kinematically connect the picture of a hadron in its rest frame with the 
picture of the same hadron in the infinite momentum frame. This is a prerequisite for any formulation 
of QCD that aims at simultaneously explaining the constituent quark model classification of hadrons 
in the particle data tables and the parton distribution functions measured in deep inelastic scattering 
processes as well as other high-energy hadronic properties that manifest themselves in collisions of fast 
moving hadrons. 



4 Scale-dependent constituent picture and the gluon condensate in hadrons 

Operators W S2Sl = U S2 W S transform a hadron state \ip) constructed in terms of quarks and gluons 
of size si into the same state but constructed in terms of quarks and gluons of size S2- Since the 
operators U s depend on interactions, they change the number of virtual particles in states they act 
on. Consequently, W S2Sl also changes the number of particles. Thus, even if one assumes that for 
si ~ s c — 1/Aqcd a proton can be represented as built from just three large constituent quarks, 
Fig. [TJ;, a change of scale from s ~ s c to s < s c is necessarily associated with creation of additional 
virtual particles whose presence is indicated in Figs. [TJd and [T^l by the large circle. Therefore, while 
the eigenvalue problem of H Sc (a Sa ) for the proton state may take the form of a Schrodinger equation 
for three constituent quarks in a potential well, the eigenvalue problem of H s (a s ) with s < s c for 
the same proton state must take the form of a Schrodinger equation for three smaller quarks that 
are accompanied by additional virtual gluons (we ignore in this discussion additional quark-anti-quark 
pairs). These gluons are condensed only in the volume of the proton. The expectation value of the 
gluon field-strength operator squared in this cloud of gluons is described below as the parameter that 
plays the role of a gluon condensate parameter in the eigenvalue problems of effective LF Hamiltonians 
derived in QCD using RGPEP. 

The title hypothesis of this article states that, although the FF vacuum is trivial, the hadronic 
gluon content is not, and it is the latter that provides the FF dynamical effects which are otherwise 
associated with the vacuum in the IF of dynamics. The remaining part of this article briefly describes 
results that support the title hypothesis. 



5 Gauge symmetry and the effective dynamics for quarks of size s < s c 

The eigenvalue problems for mesons (M) and baryons (B) expressed in terms of quarks and gluons 
from Fig. [T}d, based on the relations W SSc |12) s= = |12G% and W ss J123) Sc = |123G) S , determine wave 
functions ^> S (12G) and ^> S (123G) in the corresponding states, 

\M) S = £^.(12G0|12G>. , \B). = ^(123G)|123G) S , (9) 

12G 123G 

where the component G represents the gluons condensed in a hadron. The operator G' that creates 
the gluon component of size of a hadron from the LF vacuum is approximated by a creation operator 
for a scalar particle. The effective eigenvalue problems for mesons and baryons, H s \h) s = E\h) s , are 
projected on the basis states |12G) S and |123G) S , respectively, in order to obtain the wave functions 
ip s (12G) and ^> S (123G). The basis states are constructed to include the color-transport factors between 
i-th quark and a geometrical center, x, of quarks in a hadron, Tj = exp—ig J 1 dx^A^, and the gluon 

field operator A is approximated using the Schwinger gauge by A ll {x) — § (x — xq) v G v ^ , where xq 
denotes the position of the gluon body created by G* in the hadron and G Miy is the gluon field strength 
operator at this point (see [T51 [IS] ). 

Currently, in the absence of precise numerical information about the true low-energy effective 
Hamiltonian for LF QCD, local gauge symmetry is used to establish the structure of H s (a s ) for s < s c 
as follows [15] . Center-of-mass motion of an eigenstate of the LF Hamiltonian separates out completely 
from the eigenvalue problem and the equation one is left with is for an operator whose eigenvalue is the 
mass squared of a hadron. The operator itself is a sum of a free mass squared of constituents, denoted 
here by M 2 , plus the interactions which we do not know, say V . But we know that in the rest frame 
of constituent quarks their invariant mass squared is the square of the sum of their energies. These 
energies can be approximated by making a non-relativistic (NR) expansion, E p = m + p 2 /(2m), since 
the RGPEP vertex form factors for s < s c prevent the interactions from accelerating constituent quarks 
to large relative speeds. Thus, one can sum the individual quark energies and square the sum, neglecting 
terms smaller than p 2 . The result is a simple quadratic expression in relative momenta of quarks. This 
result is easy to obtain for arbitrary masses of individual quarks, so that one knows the right coefficients 
with which the NR momenta squared enter the invariant mass squared. By comparison of this result 
with the exact expression for the invariant mass squared in the FF of dynamics, which is quadratic in 
the transverse relative momenta but a more complicated function of longitudinal momentum fractions 
carried by quarks, one learns how to write the latter in terms of the three-momenta known in the 
NR theory using LF variables [151 . Now the crux is that we also know how gauge symmetry dictates 
interaction through the minimal coupling in the NR theory: p > p — gA. So, by analogy, we also 
know how to estimate the effects of the minimal coupling in the LF mass squared. Namely, we trace 
the consequences of the minimal coupling in the NR theory and introduce the corresponding effects in 
the LF theory. As a result, we obtain a free invariant mass squared of quarks plus interaction terms 
induced by expectation values of the gluon field in the gluon component G. Our result is the gauge 
symmetry candidate for effective H s (a s ). 

In short, the NR reasoning proceeds along the lines of Ref. [TS], except that on the basis of RGPEP 
on the LF one considers expectation values of A 2 in the gluonic component of a hadron, |G) = G^|0), 
rather than in the omnipresent vacuum |J2). After inclusion of the color transport factors, in a crude 
Abelian mean-field approximation, the expectation values of the type (G\g 2 A 2 \G) , with A = Bx r/2 
and B being the magnetic part of G^ u , render harmonic potentials of the form <PqY 2 where ip^ = 
(G|(a s /-7r)G M1/c G° [/ |G) / (G|G) and r denotes the relative position of the quarks (only their relative 
position appears because of gauge symmetry). 

By definition, in the momentum representation, the relative distances are defined as gradients with 
respect to the relative momenta. But we already know how to identify these momenta through the 
NR approximation that holds in the smallest mass eigenstates of RGPEP Hamiltonians with large s. 
Namely, in the case of mesons built from a quark of momentum p\ , anti-quark of momentum p2 (both 
having the same constituent quark mass m) and the glue component G, using notation P = pi + P2, 
x = p\ /P + , Pi = xP^ + k^, p^ = (1 — xjP 1 - — k^, one has [TS] 
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and the effective quark dynamics in the gluon condensate inside a meson takes the form 

2" 



M\ q = 4m 2 



1 



2 (KVG\ 2 1 

) 2 



2 V 3m / 2 I 9k 







For baryons one obtains (see |15j for definitions of momenta K and Q) 



M% = 9m 2 



6K' 



Q 2 +3m 



, i mpG\ 2 5 
3m / 8 



d_ 



<9Q 



(11) 



(12) 



The corresponding harmonic oscillator frequencies cum = 7t^g/(3?ti) and lub = v 5/8 u>m match phe- 
nomenologically accurate constituent quark models (e.g., see [2U] [H] [22] ) provided that one numerically 
identifies tpg with the vacuum gluon condensate parameter (fi\{a s /Tr)G fl ' vc G c lJll ,\Q) / '(f2\Q) fitted to the 
spectrum of hadron masses in the QCD sum rules [TO] E]. We reinterpret the vacuum condensate 
parameter as corresponding to the gluons condensing only inside hadrons. The ground-state eigenfunc- 
tions of operators in Eqs. ( 11 ) and ( 12 ) are exponentials of the invariant mass squared of the effective 
quarks, 



ip n = N exp 



1 




(nmy 



2nmu! 

where n — 2 for mesons and n — 3 for baryons, with cj — com and w = wg, respectively. 



(13) 



6 Observables 

Since the quarks are accompanied by the gluons, i.e., they move with respect to the gluon body G, the 
electroweak hadron form factors are obtained in the form of the Fermi motion-smeared form factors for 
the effective quarks alone. However, the smearing is a small effect, so that the resulting form factors 
closely resemble results obtained in quark models. At the same time, RGPEP appears to provide the 
right tools for understanding a transition between the asymptotic counting-rules for hard scattering 
processes and the soft, non-perturbative effects in low-momentum transfer processes. Similar comments 
apply in the case of structure functions which depend on the Bjorken x and momentum transfer Q. In 
the RGPEP, the transformation W SQSe provides a connection between the current quarks of a small 
size sq, i.e., the ones that are capable of a sudden absorption or emission of a hard-photon, and the 
effective constituent quarks of size s c ~ 1/Aqcd- This transformation is expected to describe the 
evolution of parton distributions, universally in Q and x since the size parameter controls dependence 
on the invariant masses that are functions of Q and x simultaneously. Therefore, one may expect the 
same RGPEP tools to help shed some light on the gluon saturation mechanism in QCD as seen in the 
inelastic, inclusive or semi-inclusive processes. 

It is also worth observing that the momentum variable k 1 - = / \/ 'x(l — x) that Brodsky and 
de Teramond discovered in their LF holographic AdS/QCD picture of hadrons [531 H2]> appears 
to match the one dictated here by gauge symmetry and the condensation of gluons inside hadrons. 
This matching warrants further study. Regarding this issue, we observe [15] the following: (1) The 
AdS 5th dimension appears to correspond to the quark size s in RGPEP; (2) G-induced oscillator 
potentials appear to correspond to the soft- wall (SW) models [25]; (3) the resulting identification 
K sw = 2towm = (2n/3)(fiG suggests that the soft wall in the SW models results from the gluon 
condensation in hadrons; (4) AdS/QCD SW phenomenology result km/kb ~ 1-15 ± 0.5 matches our 

result for the same quantity, (8/5) 1 / 4 ~ 1.125. These results support the idea that the condensates 
associated with the vacuum state in the IF of dynamics are actually associated only with the hadronic 
interior [57] [2S] in the FF of Hamiltonian dynamics. 

The last comment we wish to make here is that the hadron content-induced harmonic oscillator 



potential between quarks at large RGPEP size s implies LF eigenvalues M 



" 2 for states where 



r is large. This quadratic behavior of M 2 implies M 
phenomenology and string picture of hadrons. 



r, which is in agreement with the Regge 
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